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Abstract 

According to the standard cosmology, near the last scattering surface, the 
photons scattered via Compton scattering are just linearly polarized and then 
the primordial circular polarization of the CMB photons is zero. In this work 
we show that CMB polarization acquires a small degree of circular polarization 
when a background magnetic field is considered or the quantum electrodynamic 
sector of standard model is extended by Lorentz-noninvariant operators as well as 
noncommutativity. The existence of circular polarization for the CMB radiation 
may be verified during future observation programs and it represents a possible 
new channel for investigating new physics effects. 



1 Introduction 

The observation of polarization in the cosmic microwave background (CMB) and its 
correlation to the temperature anisotropics is a very promising tool to test the physics 
of early universe. The CMB radiation is expected to be linearly polarized of the order 
of 10 %. This Linearity is a result of the anisotropic Compton scattering around the 
epoch of recombination which has been widely discussed in the literature [1, 2, 3]. A 
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linearly polarized radiation is described by non-zero values for the Stokes parameters Q 
and/or U. Meanwhile the possibility of the generation of circular polarization can be 
determined by the Stokes parameter V . This parameter is usually given to be zero in the 
literature because there is no physical mechanism to generate a Stokes-V contribution at 
the last scattering surface. However the circular polarization can be generated naturally 
if a background magnetic field be present. Recently, Giovannini has shown that if the 
CMB photons are scattered via electrons in the presence of a magnetic field, a non- 
vanishing V mode can be produced [4, 5]. Furthermore, Cooray, Melchiorri and Silk 
have discussed that the CMB radiation observed today is not exactly the same as the 
field last scattered [6]. The CMB polarization may change during the propagation 
from the last scattering surface due to the presence of relativistic magnetized plasma in 
galaxy clusters. On the basis of the mechanism described in [6], the linear polarization 
of the CMB can be converted to the circular polarization under the formalism of the 
generalized Faraday rotation [7] known as the Faraday conversion [8]. The evolution of 
the Stokes parameter V given by this mechanism is obtained as 

V = 2V^, (!) 

where A(fi FC is the Faraday conversion phase shift [6]. Assuming reasonable parameters 
for the relativistic magnetized plasma in the galaxy clusters and for B = 10/xG, one can 
estimate the Faraday conversion phase shift as A<ppc ~ f ew x 10 -3 . Then with a CMB 
polarization of the order of 10~ 6 which propagates through the galaxy clusters, the out- 
going beam should contain a circular polarization of order 1CT 9 at scales corresponding 
to the galaxy clusters [6]. 

It is also possible that such an effect arise due to new physics effects. For in- 
stance another mechanism for generating the circular polarization has been reported by 
Alexander, Ochoa and Kosowsky in [9]. In this method the photon sector of quantum 
electrodynamics is extended so that photon is coupled to an external constant 4-vector 
T v via a Chern-Simons like term as 

£' — C-Maxwell + £-T 

= -\F, U F^ + ge^A^F^ (2) 

where Ct is a Lorentz symmetry violation CPT odd term and g is the coupling constant 
of interaction. Then the collision term on the right hand side of the Boltzmann equation 
is corrected by the interaction (2). The generalized Boltzmann equation considered 
in this reference is given by the quantum-mechanical evolution of the photon density 
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matrix. In fact it is shown that the time derivative of the polarization brightness 
associated with the Stokes parameter V receives a source term such that the V mode 
becomes nonzero. 

Also in reference [10], Finellia and Galavernid have considered an axion-like cosmo- 
logical pseudoscalar field acting as dark matter coupled to the photons as 

C = - 9 -±4,F^\ (3) 

where is the coupling constant between axion field and the electromagnetic 
field strength F M „ and F pu = e^ upa F pa . Then it has been shown that such an interaction 
between the pseudoscalar field and photons rotates the plane of linear polarization and 
generates the circular polarization for the cosmic microwave background. A similar 
approach with the pseudoscalar photon mixing shows a CMB circular polarization up 
to the order of KT 7 [11]. 

In this work we address new possibilities and propose new channels to produce the 
circular polarization for the CMB in three different contexts. First we consider the QED 
part of an effective field theory for Lorentz violation (LV) symmetry [12, 13, 14, 15, 16, 
17, 18] that is called the Standard-Model Extension (SME) [13]. In this case we only 
investigate that Lorentz violation terms of SME that have not been considered in [9]. 
Then we examine the noncommutative QED (NCQED) with Seiberg-Witten expansion 
of fields [19] as well as the presence of the primordial magnetic field background in 
the last scattering surface [20]. For this purpose we use the generalized Boltzmann 
equation formulated in [9], to study the evolution of Stokes parameter V under the 
influence of these new phenomena. This equation shows that a non-zero background 
field can generate the circular polarization through either the interaction of photon with 
the background itself or the corrections on the Compton scattering of photon on the 
electron in the presence of the background field. One should note that the quantum 
corrections on the scattering cross section, for the low energy particles, are usually 
negligible. Nevertheless, the circular polarization of CMB itself as a small quantity can 
be produced potentially via accumulation of these small effects. 

On the dimensional ground one can expect A<fi FC at the lowest order depends on the 
appropriate combination of the parameters of the model which is dimensionless. There- 
fore, for the magnetic field A(f) FC should be proportional to eB/m 2 that is about 10~ 19 
for the micro gauss magnetic field. Meanwhile in the NC space it should be depended 
on am 2 9 which is of the order 10~ 17 for l/y/9 about 10 TeV. For the Lorentz violation 
case the Faraday conversion phase should be linearly proportional to the dimensionless 
LV-parameters which are less than the present existing bound about 10~ 15 . Therefore 
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at first sight it seems that the effect of the magnetic field is at least two order of magni- 
tudes less than the other new interactions. But as we will show soon this is not the case 
and due to the wavelength dependence of the solutions the background magnetic field 
has main contribution on the polarization of the CMB. However, it should be also noted 
that the resulting circular polarization modes are certainly small due to the smallness 
of quantum effects. Meanwhile the probe of circular polarization of the CMB radiation 
through the future subtle experiments, certainly give more information about the new 
physics or the primordial magnetic fields in the early universe. The main condition 
which is assumed everywhere during the paper is a cold plasma at the last scattering 
surface at a red shift of about 1100 or an age of about 180,000 (f2 /i 2 ) -1 ^ 2 yrs and a 
temperature of the order of 0(eV). The low temperature is always applicable for the 
calculation of the Compton scattering amplitude. 

The paper is organized as follows: In section 2 we review the Stokes parameters 
formalism and the generalized Boltzmann equation. In section 3 the generation of 
circular polarization due to the Lorentz violation is discussed. In section 4 the same 
calculation is done for the noncommutativity. In section 5 we examine the generation 
of circular polarization in the presence of a background magnetic field. Finally in the 
last section we discuss about the result. 



2 Stokes parameters and Boltzmann equation 

The polarization of CMB is usually characterized by means of the Stokes parameters 
of radiation: I, Q, U and V [3]. For a quasi-monochromatic wave propagating in the z- 
direction, in which the electric and magnetic fields vibrate on the x-y plane, the electric 
field E in a given point can be written as 

E x = a x (t) cos(u t - S x ), E y = a y (t) cos(w t - S y ), (4) 

where the wave is nearly monochromatic with frequency uj . Then the Stokes parameters 
are defined by time averaging of the parameters of electric fields (4) as follows 

(5) 
(6) 

Sy)) , (7) 

<y> • (8) 

These parameters are physically interpreted as follows: I is the total intensity of wave, Q 
measures the difference between x and y polarization, the parameter U gives the phase 



I = (al) + (a 2 y ) , 
Q = (a 2 x ) - (a 2 y ) , 
U = (a x a y cos(5 x - 
V ee (a x a y sin(4 - 
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information for the two linear polarization and V determines the difference between pos- 
itive and negative circular polarization. The time evolution of these Stokes parameters 
is given through Boltzmann equation. We now turn to the Boltzmann equation which 
is vital in our study. The CMB radiation is generally described by the phase space 
distribution function / for each polarization state. Boltzmann equation is a systematic 
mechanism in order to describe the evolution of the distribution function under gravity 
and collisions. The classical Boltzmann equation generally is written as [21] 

f t =C\fl (9) 

where the left hand side is known as the Liouville term deals with the effects of grav- 
itational perturbations about the homogeneous cosmology. An important point that 
must be mentioned here is that we have disregarded the scalar fluctuation of metric on 
the left hand side of the Boltzmann equation. The right hand side of the Boltzmann 
equation contains all possible collision terms. For photons, the Compton scattering with 
free electrons, j(p) + e(q) -B- j(p') + e(q'), is more important and must be included on 
the right hand side. To compute the evolution of polarization we follow the approach 
of references [3, 9]. At first, the distribution function / is generalized to the density 
matrix 

P = 

= ^(It + Qa 3 + Ua l + Va 2 ), (10) 

where 1 is the identity matrix and <7j are the Pauli spin matrices. The density matrix 
Pij is related to the number operator A?(k) = aj(k)aj(k) as 

<P y (k)> = (27r) 3 2fc°^ 3 )(0K-(k). (11) 
and the time evolution of the number operator is given by 




;Aj) (*) - i([#int(*), Aj]> - I dt'([H in t(t - t% [H iQt (t')Aj}})- (12) 



d 

dt' 

Hence in terms of p, the evolution equation leads to 



DC 



(27r) 3 5 3 (0)2A; ^(0,k)=z([^ int (0),X) l ,(k)])-i / ,//' ([H iM (t') [HiM- A y (k)]]) . 
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(13) 

where the first term on the right hand side is referred as the refractive term and the 
second one as the damping term. Note that in (13), all factors of 5 3 (0) will be canceled 
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from the final expressions. Equation (13) can be viewed as a quantum mechanical 
Boltzmann equation for the phase space function p. In the case where the interaction 
term is the convenient QED interaction, the time evolution of the Stokes parameter V 
is always equal to zero which leads to the absence of circular degrees of freedom for the 
CMB photons. In this paper during the next sections, we will show that if one consider 
beyond standard model interaction terms, the time evolution of the V receives a source 
term which is interpreted as a non-zero circular polarization. If the Stokes parameter 
V be non-zero, then it contributes on the new angular power spectrum elements such 
as Cf T , Cj E and Cj v where may be detectable in future CMB experiments. 

3 Generation of circular polarization in the pres- 
ence of Lorentz violation terms 

In the present section we consider a new class of Lorentz invariance violation (LIV) 
terms as a generic class of interaction between photons and an external field and ex- 
plore whether they can produce circular polarization for the CMB photons or not. 
Although the source of these asymmetry terms is unknown, physicists usually suppose 
that they come from presumably short-distance physics which manifest themselves in 
the interactions of standard model (SM) fields [12, 13, 14, 15, 16, 17, 18]. The renor- 
malizable sector of the general QED extension for a single Dirac field ip of mass m and 
a photon also has been studied by some authors for which the Lagrangian C is written 
as 

C = £ QED + CtTv ron + £f/y° n , (14) 

where 

£qed = i h*Vpi> - ™W - \f^\ (15) 

electron = « C ^ V ^ + » d^T^, (16) 

photon = _ 1 {kF )^ F ^ F ^ + ^(k AF ) a e a ^F^, (17) 

where the X> M is the covariant derivative. In the electron sector, c MJ , and are the di- 
mensionless Hermitian coefficients with symmetric and antisymmetric space-time com- 
ponents. For the photon sector, it will be useful to decompose the coefficient k F into 
a tensor with 10 independent components analogous to the Weyl tensor in general rel- 
ativity and one with 9 components analogous to the trace-free Ricci tensor with the 
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following symmetries of the Riemann tensor 

(kFJfivaP — — (kFJvfiaP = ~ (k F )^ u pa = (k F ) vfM /3 a = (k F ) a p^u- (18) 

So it contains 19 independent real components. Various experiments have bounded 
these coefficients both in fermion and photon sectors. For instance, the only rotation- 
invariant component of {kp) l i V pa is constrained to < 1CT 23 [12] and all other components 
of (kF)fMu/3a which are associated with violations of rotation invariance can be bounded 
of about < 1CT 27 [13]. The high-quality spectropolarimetery of distant galaxies at 
infrared, optical, and ultraviolet frequencies, also constrains the coefficients hp to less 
than 3 x 10~ 32 [15]. In references [16, 17], another bound has been estimated. For the 
couplings c M „ and d^ u in the electron sector, one can find a variety of bounds such as 
Cfj,u < 10~ 15 and d^ v < 10~ 14 in [18] and references there in. 

The \{kAF) a tafiiivA^ term in the photon sector, has been studied in [9] and it 
was shown that the CMB photons acquire small degree of circular polarization. In this 
section we consider the three remaining terms in the LIV Lagrangian. We study the 
contributions of each terms in the Boltzmann equation and calculate the evolution of 
the Stokes parameters. Firstly, the photon sector term with hp coupling and then the 
electron sector with c M „ and d^ v couplings are investigated. 



3.1 The pure photon sector 

In this part the Lorentz violation terms (17) are investigated. These terms do not 
modify the Compton scattering but the dynamics of the CMB photons can be influ- 
enced. Here we consider H int = — \ {kp) ilva ^F OL ^F ilv on the right hand side of the 
generalized Boltzmann equation (13). Using the canonical commutation relations of 
the creation and annihilation operators and their expectation values given in [3], the 
equations (4.11a)-(4.11e), one can calculate the refractive term as 

n int (t), Ptf(k)] } = ~ (k F r^(2n) 3 5^(0)(k a k,e* s , l3 e su ) [<Wv,(k) - 5 jV p is (k)] , (19) 

where e sll (k) are the photon polarization 4-vectors with s = 1,2. Now using the sym- 
metry properties of k F va ^ summarized in Eq. (18), one can write as 

j t p tJ (k) = -1 (^(U^) [W-i(k) - <W«(k)] • (20) 

The photon density matrix can be expanded about a uniform unpolarized distribution 
as follows 

PV = P$ ) +P§\ (21) 
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where p^ = p 22 and pf 2 = P21 = 0- Therefore the components of the time derivative 
of the density matrix are given by 

p$(k) = ~ {k^Kk^elJS ~ ^e\JSl (22) 
pS(k) = ~ {k F r^k a k,[e w e\JS - ei&JSl (23) 
p«(k) = --1 (k F r^k a k,[e w e* lu p^ + e ip e* 2 J£ - E ip e* 2 J$ - e 2p e* 2 J$l (24) 

p%(k) = ~ {k F Y^k a k,[e 2 pe\JS + ew&p® - e^e\JS - e^e\J 2 X 2 \ (25) 

Note that the polarization vectors e s ^(q) are real. Hence the evaluation of the Stokes 
parameters are given as 

J (1) = 0, (26) 

1 fi 

Q (V) = "To (k F r^k a k,(e w e* 2u - e w e\ v ) (27) 



k° 
4 



= - A (M^aM^L - e ip e* 2v ) Q (1) , (28) 



= (M^^{(^L - - + e^L) (29) 

Therefore in the first order of Lorentz violation parameter k F , V is nonzero. This 
shows that although the usual Compton scattering of CMB photons can't generate 
circular polarization, the interaction of photon with a nontrivial LIV background clearly 
produces a nonzero circular polarization. To have an estimate of the Faraday phase 
conversion we rewrite (29) as follows 

VW = ik F T^{L% v [/« - L% v Q«}, (30) 

where L u a ^ v = ff |fK e i/3 e L - e 2^* 2u ) and a corresponding relation for L^ v . Now for 
Ao = 1cm, z=1000 and for 1 kpc one has 

A(f) FC ~ 4 x 10 20 A; F , (31) 
or to have a phase of order one we should have k F ~ 1CT 20 . 
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3.2 The electron sector 



3.2.1 & v term 

We now turn to the LIV terms in the electron sector (16). Calculating the refractive 
term on the right hand side of the generalized Boltzmann equation (13), leads to V = 
for interaction (16). Therefore, We move on to the damping term. This term gives 
corrections to the Compton scattering of the CMB photons and electrons. These cor- 
rections modify the right hand side of the equation (13) and then may influence the 
evolution of the Stokes parameter V. In this part we investigate whether these correc- 
tions lead to a circular polarization mode or not. Firstly, the gauge field is expanded 
in terms of annihilation and creation operators a and and the fermion fields ip and 
■0 are expended in terms of b and b\ Then the interaction Hamiltonian can be written 
in the momentum space as 

H int (t) = J dqd^dpdp' (27r) 3 5( 3) (q / + p / -q-p)exp(2t(g' +p' -g°-p )) 

x b\,{q')a\ l {p'){M l + M 2 )a s {p)b r {q), (32) 

where q and q' are the momenta of electrons and p and p' are the momenta of photons. 
Also we have used the following shorthand notations 

d 3 q m d 3 p 
rfq= (2^V rfp= (2^2A (33) 

The Mi and M. 2 are the amplitude for the Compton scattering process up to the first 
order of LIV parameters d J,u and d^ v . Consequently, the commutator in the damping 
term, finds the following form 

H hlt (t),V tJ (k)^ = J dqdq'dpdp' (27r)V 3 )(q / + p / -q-p)(A^ 1 + A^ 2 ) 

x[bl l (q%,(q)al l ( P y s (p)2p (27r) 3 5 ts 5 :i (p - k) 

-^(g')^(g)4(p / )a s (p)2p /0 (27r) 3 5, s ,5 3 (p / - k)]. (34) 

In appendix A we detail the calculation of (34) and it's expectation value. Quoting 
the results of Appendix A, and expanding the density matrix as = pf) + p{^\ the 
evolution of Stokes parameters is derived as 

J (1) = 0, (35) 
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Q (1) = ^o/ rfqne ( q ) ~J { q - ke y» ~ k ^ e \ ~ <l- k€ l € l + Kq-e^l 



Ac" u 



q v (q-z% + q-e 2 el)} + J^^% + Kk^q.^q.e 2 ) J V W , (36) 



uw = 



2k° 



rn 



dqn e (q) 



2& v 
q.k 



{QM e ^u ~ e M - k „q-e S + Kq-e 



2c txv 1 
+2q u (q.e 1 el - g.e 2 ^)} + ——{q^ + Kk^q.^q.e 1 - q.e 2 q.e 2 )\ V (1) , (37) 



(q.k) 



and 



le 



«2c fMU 
—j^{(q-k( e Wu + 44) - k v q.ee\ + Kq.t 1 ^ 



"2& v 

- 1 - 2 2 11 ' ;(q,q, + Kk,)(q.e\e 2 ))(-Q^ 



+qv{q-e 1 el + q.e 2 el)} + 



(q.kf 



+ {^{q-H e l e l - 44) - k v q.t l t\ + k v q.e 2 t\ 



+g„(g.e 1 4 - g.e 2 ej)} - 7^(<7^ + M^Xg.eVe 1 " g.e 2 g.e 2 ) ) f/« }. (38) 

An exact solution for the evolution of the Stokes parameters involving the coupled 
evolution equation for all of the species can only be accomplished through numerical 
integration. The time derivative of the Stokes parameter V has a source term which 
indicates the appearance of circular polarization for the CMB photons if they scatter 
from electrons through Compton process with an extra Lorentz violation c^ u term. Here 
again to have an estimate of the value of circular polarization we can find the order 
of magnitude of the Faraday phase conversion. For this purpose we consider only the 
largest term and confine ourselves to q ~ m 3> |q|, |k|. Therefore one has 



VW = 



ie 2 m _ 
-2kJ V ' 



(? ,/ Vi l v l/ {-v.6 1 v.6 2 Q( 1) + 2{v.e 2 v.e 2 - u.eW)l/ (1) }, 



where 



f_f_q 
J (2vr 



n P 



f d 3 q 

; J j7^qiM<i) = mviUe, 



(39) 



(40) 



and v = (l,v). Thus for A = 1cm, z=1000 and for length scale about 1 kpc and 
number density of electrons about 0.1 per cm 3 one has 



A0 FC ~ 1 x 10 2 Vl, 



(41) 
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or to have a phase of order one we should have \d iv \ ~ 1CT 20 . One should note that 
if we take the intergalactic electron number density of about 1CT 8 per cm 3 then for a 
length scale about few Gpc the result only changes an order of magnitude. 

3.2.2 df™ term 

We now ignore the and move to the d^ u term of the Lagrangian (16), 

£ = % -^V^ - + l - dT^lP^- (42) 

Similar to the case, one can evaluate the Compton amplitudes Aii and hA 2 and 
then the interaction Hamiltonian in order to find the damping term. The computations 
are detailed in appendix B. Using that results, the evolution of Stokes parameter V is 
given by 

V {1) = 0. (43) 

Hence the d^ v correction to the electron sector does not generate any circular component 
for the polarization of CMB. Although considering the higher order corrections of d^ u 
may result in a non-zero source term for the V. 

4 Circular polarization due to noncommutativity 

In this section we consider the noncommutative corrections to the QED interaction term. 
There is two fundamentally different approaches in which one can find a closure form for 
the Lie algebra of a noncommutative gauge theory [22, 23, 24, 25, 26, 27, 28, 29, 30, 31]. 
In the approach proposed by Chaichian et al. [22], the group closure property is given 
only for U(N) gauge theories with matter content in the fundamental or adjoint rep- 
resentations. This kind of model building has been criticized because it leads to some 
conclusions for the standard model which is not consistent with the our knowledge 
of particle physics. Wess and his collaboration have developed another approach to 
construct a consistent noncommutative gauge theory [23]. In this approach, the Lie 
algebra is extended to an enveloping algebra. The mapping between the noncommu- 
tative gauge theory and the effective theory on the usual commutative space-time is 
obtained by requiring that the theory be invariant under both noncommutative and the 
usual commutative gauge transformations. These requirements are given through the 
solutions of some differential equations known as the Seiberg- Witten map (SWM) [23] 
appeared originally in the context of string theory. The basic property of the SWM 
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Figure 1: The fermion-photon-fermion vertex of the noncommutative 
QED: ieQf {7^ - \ [{Pout0p in )l^ - (PoutO)^in- m f ) - (fe^^r ™/)]} 



method is that the new interactions are added but no additional particles are intro- 
duced. Also in this approach both the gauge field and the gauge transformation 
parameter A are expanded in terms of noncommutative parameter 

A, = A, + \e^{d a A, + F a „ Ap} + 0(9 2 ), 

A = A + ^{d a A,Ap} + 0(9 2 ). (44) 

The advantage is that this construction can be applied to any gauge theory with arbi- 
trary matter representation. The SWM formulation of the noncommutative standard 
model thus uses the standard model gauge group. Following [23], the noncommutative 
generalization of QED sector is defined by the action 



d x 



tp *i(p* rji) - mi) ★ i) - -IrF^F^ 

2 



(45) 



where 

F^v = + ^{F m , F u/3 } - ^{A a , {Op + V P )F IW } + 0(9 2 ), (46) 

and 

$ = 4 + X -e9^ A^ + 0(9 2 ). (47) 

The Feynman rules for the fermion-photon-fermion vertex is modified as what given in 
Fig. 1. The noncommutative parameter 9 in SWM approach has been bounded through 
different experiments. Using the experiment data, the best bounds are of the orders of 
1 - 10 TeV [32]. 

Now we turn to the Boltzmann equation and perform the damping term. Using the 
new vertex, the amplitude of the Compton scattering on a noncommutative background 
is modified and gives corrections. Details of the calculation of the evolution of density 
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matrix elements are presented in Appendix C. Quoting the results (94-95) and similar 
to the previous section, the evolution of the Stokes parameter V is given as 



%(? f ( 

= 7^7-^ J dqn e (q)\2(q9e 1 q • e 1 + q9e 2 q ■ e 2 )(p 12 - P21) ~ 2(q9e 2 q ■ e x + q9e x q ■ e 2 )p u 

+2(q9e 2 q ■ e x + qde iq ■ e 2 )p 22 ) , (48) 
where it can be more simplified as follows 

• 2 r 

V {1) = ]T^J dqn e (q)(q9e 2 q ■ e x + q9e x q ■ e 2 )Q (1 \ (49) 

Hence similar to the Lorentz symmetry violating interaction terms discussed in the 
previous section, the noncommutative quantum electrodynamics based on the SWM, 
also can generate the circular polarization for the CMB radiation in the first order of 
noncommutative parameter 9. Now (49) can be cast into 

V {1) = % 4^(v.6.e 2 v ■ e 1 + v.e.e.v ■ e 2 )Q«, (50) 

where A^c is the value of the non-commutativity parameter and 9 is unit vector that 
shows the non-commutativity direction and is a constant. Therefore the Faraday phase 
conversion for A NC ~ 1 — 10 TeV will be around 10~ 7 — 10~ 9 radians for A = 1cm, 
z=1000 and for length scale about 1 kpc. 



5 Circular polarization and magnetic field 

When photons impinge electrons in the presence of a magnetic field, the circular po- 
larization is generated naturally which leads to a non vanishing Stokes parameter V 
whose magnitude depends upon the intensity of the radiation field [4, 5]. For instance, 
If the magnetic field is considered of the order of B ~ nG, the angular power spectra 
for small values of I are derived as £>(10- 7 )(/iK) 2 and £>(10- 15 )(/iK) 2 for the VT and 
VV components respectively [4, 5]. In these references, the electron-photon scattering is 
computed by classical considerations. Here we want to investigate a similar phenomenon 
by considering the Compton scattering of electrons and photons. A physical situation is 
considered when prior to recombination a cold plasma of electrons is supplemented by a 
weak magnetic field whose typical inhomogeneity scale is at least comparable to Hubble 
radius H~ l at the corresponding epoch. In the epoch of recombination, magnetic field 
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is estimated to be about 10~ 3 /iG to 10 _1 /xG [33]. This value of magnetic field is very 
weak compared to the critical value B c = m 2 e /e = 4.414 x 10 13 G. For this reason during 
our strategy, the effect of magnetic field on the electron wave functions is ignored and 
only the electron propagators are corrected. The propagator of a Dirac particle in an 
external constant electromagnetic field background is well known from Schwinger work 
[34] 

S(x, y) = y) J ^ e HP.(x-y)5 (p)> (51) 
where the phase factor $(x,y) is defined by 

^(x,y) = e X p(- l -^x a B a py^j, (52) 

in which B a = — jB^x 13 and B a p represents the strength background magnetic field 
tensor. The other part of Schwinger propagator, S(p), can be expanded up to the linear 
background field as 

S(p) = S (p) + S B (p) + O(B 2 ), (53) 

for which So is the usual electron propagator and Sb is the linear modified part of the 
Schwinger propagator given by [35] 

1 l a (j> — m)^P 



S B {p) = --eB aP >J J> . (54) 



Therefore invoking the above considerations, the Compton scattering matrix element 
in the presence of background magnetic field is given by 

M = Mi+M 2 , (55) 

where 

M 1 = u r >(<f) J d 4 xd 4 x'e ix '- q '(-ie£ s ,(p')e lx '- p ')S(x',x) (-ie{ s (p)e~ lx - p ) e- ix - q u r (q), (56) 
and 

M 2 = u r ,{q') J d 4 xd 4 x'e ix '- q ' (-ie/t s (p)e- ix '- p ^ S(x',x) {-ie^{p')e lx - p '^j e""-% r (g), 

(57) 

where u r (q) and u r >{q') respectively are the wave functions for the incoming electron 
with spin index r and outgoing electron with spin index r' . Also e^{p) and e^(p') 
are polarization vectors for incoming and outgoing photon chosen to be real and with 
indexes s,s' — 1,2 labeling the physical transverse polarizations of photon. Using the 
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linear expanded Schwinger propagator (53) and by expanding the $(x, y) factor up to 
order B, for the amplitude (56) one can write 

M 1 = u r ,(q') J d 4 xd A x'-0^e ix '^' +p '-^e-^ +p - k ^-ie^i P '^ 

x (s (k) + S B (k) - jx' a B a ^So(k)^ (-ie/t s (p))u r (q). (58) 

Integrating over the first two terms of (58) gives 

{2^5\q'+ V '-q- V )uAq') {-ie^ip')) (S (q + p) + S B (q + p)) (-ie/t s (p)) u r (q) . (59) 

The remaining |i ,a B tt ^i^o(fc) term explicitly violates the four momentum conser- 
vation. One can assume that during the interaction time a small momentum Sq^ is 
transferred from background into the internal electron of the associated Feynman di- 
agram so that the momentum conservation will be preserved and takes the following 
form 

q» + p'» = q» + p» + (60) 

Substituting the condition (60) into the last term of (58) and expand the exponential 
factor yields 

l -^B aP j d 4 xd A x'e i{x '- x) - iq+p ' k) x' a x p {l + ix'.Sq). (61) 

In the above expression, the term that is independent of 5q is vanished because of anti- 
symmetric properties of B a/ 3 with respect to a and (3 while the term that is dependent 
on 5q remains. From classical electrodynamics we know that the four-momentum of 
electron that propagates between two space-time points, changes similar to the Lorentz 
force Apu = —eB llv Ax v . Then the magnitude of Sq^ can be approximated to be of the 
order of external background field 5q ~ eB. Therefore this term becomes of the order 
of 0(B 2 ) and hence can be ignored in the scattering amplitude With a similar 

procedure, the amplitude M. 2 is given by 

M 2 = (27r) 4 5 4 (q'+p'-q-p)u rl (q r ) (-iejk 9 (p)) (S (q - p') + S B (q- p')) (-iefatf)) u r {q). 

(62) 

The ^-independent terms of M\ and M 2 give the usual Compton scattering. The 
correction up to the linear background field is represented as M. B and is written as 

M B = Mib + M 2B 

= (2nf5\q' + p' - q -p)u{q')[{-ie^{p')) S B (q + p) (-u&(p)) 

+ (-iejt s (p)) S B (q - p') (-ie^(p')) 1 Mq)- (63) 
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Using these amplitudes, the damping term on the right hand side of the Boltzmann 
equation can be computed. Computing the damping term is detailed in appendix D. 
By using these results the density matrix element and then the evolution of the Stokes 
parameter V up to order of e 4 is derived as 

V (1) =t(p 12 -p 2 l) 

xn e (q) [((qMp)) 2 + (g.e 2 (p)) 2 )(/ (1) (k) - / (1) (p)) " (^i(p)^i(p) ~ qMp)qMp))Q 
-2q.e l {p)q.e 2 {p)U {l \v)} +0(k,p), 

where = —eB^ u q u . Hence similar to the LIV interaction terms and the noncom- 
mutativity, corrections due to the external magnetic field can generate the circular 
polarization for the CMB radiation as well. But in contrast, it is not needed a prior 
linear polarization to have a nonzero Stokes-V parameter that is in agreement with 
references [4, 5]. But it should be noted that in [4, 5] the scattering matrix is classi- 
cally obtained in the electron-ion magnetized plasma, while in our model the Compton 
scattering is calculated by considering the correction on the electron propagator in an 
external magnetic field. Here one can easily see that the circular polarization is linearly 
proportional to the background magnetic field. Meanwhile if we integrate (64) over the 
momentum p and use the approximation q ~ m ^> |q|, |k| then one can see that the 
correction will be proportional to A 3 as follows 

•m e 4 mA 3 f dVt . ~ A ,,,,eB, „ r TT _ B uu , 

V (1) = ^— / ^n e (v.(k - p)) 2 ( - F t;, p, e u e 2 , J, U, Q, , (65) 

57T J 47T m z B 

where F can be easily define by comparing (64) and (65). Here again one can estimate 
the Faraday phase conversion for z = 1000 as 

A^ lrarf (_£_) (J_). (66 ) 

6 Discussion and Conclusion 

In this work we have considered the evolution of the Stokes parameters in the context of a 
general extension of quantum electrodynamics allowing for Lorentz violation (LVQED), 
NCQED and in the presence of the background magnetic field. In fact all nonzero 
backgrounds can potentially produce the circular polarization for the CMB radiation. 
In the LVQED we have explored three cases, one case in the photon sector and two cases 
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in the electron sector. In the former one the induced circular polarization is provided 
through the interaction of photon with the background Lorentz violation field, Eq.(30), 
while in the latter cases the corrections on the Compton scattering in the LVQED 
cause the polarization, see Eqs.(39) and (43). It is shown that the Stokes parameter 

V in both sectors are linearly proportional to the dimensionless Lorentz parameters &f 
and c that is in contrast with the result of [9] where the second term of Lagrangian 
(17) is considered as the interaction term and it is found that V depends quadratically 
on the Icaf that is too small with respect to our result. We have also estimated the 
Faraday phase conversion in the LV case to give an upper bound on the dimensionless 
Lorentz parameters c and kp about 10" 20 that is comparable to the current bounds on 
these parameters. In the NCQED we have found the nonzero parameter V through a 
correction on the Compton scattering in the NC space that is linearly proportional with 
the parameter of non commutativity as given in (50). Our estimation on the Faraday 
phase conversion is about 1CT 7 — 10~ 9 radians for A^c ~ 1 — lOTeV that is too small 
to be detectable. Finally we have examined the Stokes parameter V for the background 
magnetic field. In this case the refractive term is zero but the damping term has a 
nonzero V parameter that depends linearly on the magnetic field, see Eq.(65). This 
result can be compared with the result of [6] in which the circular polarization of CMB 
due to the Faraday conversion depends quadratically on the background magnetic field. 
Furthermore, in contrast with [6], all the other Stokes parameters (I, Q and U) have some 
contributions on the generation of the Stokes-V parameter in the background magnetic 
field. In fact it is not necessary a prior linear polarization to have a nonzero Stokes-V 
parameter that is in agreement with references [4, 5]. By comparing our results together 
we see a A 3 -dependence for the magnetic field and LV on one hand and A-dependence 
for non-commutativity on the other hand which leads to different spectrum for different 
wavelengths. Because the required sensitivity level to detect circular polarization of 
CMB is beyond the current instrumental techniques, the experimental bound on the 

V modes is rare. The bound reported in [36], constrains it as V < 1 /iK. However, 
the post-Planck CMB polarization experiments are already under study [37]. Then it 
is time to address the question whether CMB is circularly polarized or not and what 
physical information can be extracted from its measurement. Also it seems that if a 
new experimental device can be built with sufficient polarimetric sensitivity beyond 
the /xK range and with good control of systematic errors, the resulting observations 
may have the potential to realize the physics at energy scales far beyond the reach of 
the largest feasible particle accelerators and provide a window to new physics at these 
energy scales. Therefore we encourage experimenters to make experiments in order to 
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measure the circular polarization of cosmic microwave radiation. 



A The calculation of density matrix elements with 
& v coefficients 

In this appendix at first the scattering amplitudes M.\ and .M2 are calculated as follows. 
For the Lagrangian (16) with & v term, the equation of motion for the electron if) is found 

as 

(ifi - m + i c v ^ v d^i\) = 0. (67) 
for which the modified Dirac equation is written as 

{4-m+ c^q^uM) = 0, (68) 

in which ip = u s (q) exp(-iq.x). In the Weyl chiral representation for the Dirac matrices, 
after some straightforward calculations, one finds the following solution for the Dirac 
equation (68) 

u.M)={ V^ + ^'M (69) 

where a = (1,ct) and a = (1, — a). The electron propagator can be also calculated 
straightforwardly with a correction in the denominator as 

S(p) = . (70) 

i - m + c^q^ u 

Since |c^ M | 1, the propagator can be expanded up to the linear order of c Ufl as 

S(q) = S (q) + S (q)(ic^q^)S (q) + ■■■ 

+ ' {lc ^ %lv) ' + ... ) (71) 



§ — m § — m § — m 

where S (q) is the usual electron propagator. The u s (q) solutions can be also expanded 
in terms of c Ufl , 

« a (5)=«o(?)+«(5) + ---, (72) 

where uo(q) = u s o(q) is the free electron solution for the usual Dirac equation in the 
absence of c Ufl term and u(q) = u s (q) is the correction of the c Ufl order. The dots show 
the higher order corrections. Consequently, the Compton scattering amplitude M. is 
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modified with corrections given by 



M = Mi + M 2 

= u(q')(-ie) (Y + c^ lu ) e; s '(p')S(q + p)(-ie) (Y + c^ lu ) e s p (p)u(q) 
+ u(q>)(-ie) (Y + c»o lv ) e s p (p)S(q - p')(-ie) (Y + c^ lv ) (p')u(q), (73) 

where e(p) and e'(p') are transverse polarization vectors for the external real photons. 
The scattering amplitudes up to the order of c Ufl are simplified as follows 

M 1 = -e 2 u (q'){Ye*; '(p')S (q + p)c^ lv e s p + Y^' (p')S (q + p) {ic^q plu ) x 
So(q + ph%(pW;'(p')S (q + ph%(p)}u (q) - 

e 2 u(q'h%(p)So(q + P )Ye s p u(q), (74) 

and 

M 2 = -e 2 u Q {q'){Ye s p {p)S Q {q - p' W;'(p')Ye s p (p')S (q - p') (u?% lv ) X 

So(q -p')Ye*;\p') + c^^e s p (p)S (q-p')Y^'(p')}u (q) - 

e 2 u(q')Ye s p (p)S (q-p'h% s, u(q). (75) 

Substituting these amplitudes in the Eq. (34) and using the expectation values computed 
in [3] and using the photon polarization vector properties €i(k) • k — and ej • €j = 5^, 
it follows that 

n mt (t),V t ^)]) = e 2 (27r) 3 <5 3 (0) J dqn e (q) [S siPs/j (k) - 5 s , jPsi (k)} 

x ([^){c^e s p {k){q.ke s ;{k) - q.e* s 'k u ) + c»% s ' {k){q.ke s v {k) - q.e s k u ) + 

c^q p {q.e s el s ' + q.e* s 'e s u )} + |^(» + k p k v ){q.e s q.e* s ')) , (76) 

where the quantity n e (q) is the number density of electrons of momentum q per unit 
volume. It is now straightforward to calculate the time derivative of the density matrix 
which result in 

e 2 f r &" v 

Pn(k) = -^—^ J dqn e (q) [ — [{{q-ke^el - k v q.t 2 e 1 ^ + {q.ke l v e 2 p - k v q.e l e 2 p ) 

+qu(q-e 1 e 2 p + g.e 2 ej)}p 21 - {(g./ceje* - k v q.e l e 2 p ) + {q.ke\e 2 v - k v q.e 2 e x ^ 

+q v {q-e 2 el + q.eh^jpu] + + k tl k v )q.e l q.e 2 (p 2 i - Pi 2 )]J , (77) 
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e f r c^ v 

+q u {q-e 2 e\ + q.eh^jpu - {{q.kelel - Avg.e 2 eJ) 
+(Q-kelel - Kq.e 1 ^) + q v {q.e x e\ + g.e 2 ej)}p 2 i] 



2c/ 4 " 



{(<? M ^ + k iJ k v )q.e l q.e 2 {p X2 - P21)} 



(78) 



Pi2(k) 



2mA; 



dqn e 



(q)^[{(?-^i e ^ ~~ M-e 2 ^) + (g./ce^ej - Kq^e^) 



+Qu(q-e + g.e e M )}(p 22 - p n ) + {g.fc(e M e„ - e^ej - fc„g.e e M + k v q.e 
- (M-e^i - k v q.e 2 el) + 2q v {q.e x e\ - g.e 2 e 2 )}pi 2 ] 

2c M " 

+ 7^p{(» + k fl k u )[q.e l q.e 2 (p 22 - p n ) + ((g.e 1 ) 2 - (g.e 2 ) 2 )p 12 ]} 



(79) 



and 



P2i(k) 



2mA; 



dqn e 



(^) — [{(Q-ke 2 A- k ^T 



+(q-kelel - Kq.e 2 el) + g M (g.e 2 e* + g.e 1 e 2 )}(pn - P22) 

+{2q.k(e 2 p el - eje*) - 2{k u q.e 2 el - Kq.t 1 ^) + 2q l/ (q.e 2 e 2 i - q.e 1 e\)} p 2 i\ 



2& v 



(?f 9" + ^k u )[q.e 1 q.e 2 (p 11 - p 22 ) + ((g.e 2 ) 2 - (g.e 2 ) 2 )p 21 ]j . (80) 

B The calculation of density matrix elements with 
coefficients 

In this appendix we compute the density matrix elements with interaction given by 
(42). The equation of motion for the Lagrangian (42) is calculated as 



u s (q) = 

and the propagator changes as 

S(q) 



y /q.a + {d^a u ) ( s 
yjq.a - (d^a u ) C 



(81) 



j) - m + d u ^p/y b ^ v 



(82) 
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similar to the &- v case, one can write the following expansion for the propagator up to 
the order 0{dT) 



S{q) = S (q) + S (p)(id^q p Yl»)S (q) 



+ (id^q^h • (83) 



m a — m a — m 



Then the Compton scattering matrix elements M\ and M.2 in this theory are calculated 

as 

M = Mi + M 2 

= u{q'){-ie) (Y + d^Ylu) e*;'(k')S(q + p){-ie) (Y + d^VS) ^(p)«(?) 
+ u(q')(-ie) (Y + d^Ylu) e s p (k)S(q - p')(-ie) (Y + d^Yl.) ' (pk')u(q) ,(84) 



where 



M^i-ffaffli-ftfWoiq + p) {d up l 5 lu) e s p 

+ Ye; s, (k')S (q + p) (id»%Yl») S (q + p)Ye S p (p) 
+d^Yl,e*;'(p')S (q + p)Ye s p (p)}u (q) 

-e 2 u(q'h%(p')So(q + p)Ye s p u{q), (85) 



and 



M 2 = {-e 2 )u Q {q'){Yt s p {p)S Q {q - p') {d^Yl») <f V) + 

Ye s p (p')S (q - p') {id^qpYlv) S (q ~ p')Y^\p') + 
d up l 5 l,e p (p)S (q -p')Ye* p s \p')}u Q {q) + 

(-e 2 )u(q')Ye s p (p)S (q - p')Y^ '«(?). (86) 

Substituting these amplitudes in the Eq. (34), and using the expectation values calcu- 
lated in [3], leads to 

o2 



2 f 

n int (t),V l3 (k)]) = -^(2tt) 3 5 3 (0) dqn e (q)[5 isPslj (k)-5 js/ p is (^)} 

J / ^ 777- j 

d v » 
Wq) 



{ (4 S ' 4 - 4 S '4)[( _2 ^^(^^ + fc /* fc ") ~~ 2 Q<xQ P (Qi*k v + k p qu)) 

1d v ^ 

K s 'q a (m p - q p el) - e^q a {q^' - q p ef)\ - —^P^g^}, (87) 
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where e a/3crp is the Levi-Civita symbol appeared through tracing process when 7 5 is 
present. Similar to the c^ u case, the time evolution of the pu and p 2 \ are obtained as 

Pi2(k) 



and 



J dqn e (q){/ 



\Aid vp 



2mk° I "^'^ l Pl2 [~qJc~ ( e a e P ~ V P )1 + (P22 ~ pn) 

2d^ 



\ q .ky 



+ (k 2 q a K - 2m 2 q^k a + 2{q.k)q a q^)e al3 u P } + —re a ^ u 
[Qaq P (elel - e\e\) + q a qp{e]/ p - ejej) + 2iq ii q$<? e 1 p \ }, 



(88) 



p 2 i(k) = 



2mA; 



dqn e (q)q°\ - p 2 i 



Aid"' 1 
q.k 



-{e\t\ - elel)[{-2q a k p {q^q v + k^k v ) - 2q a q p {q p k v + k^qu))e al3pa 



2d Uf * 



+ (k qjip ~ 2m q^k a + 2(q.k)q a q fl )e ap „ } H r e ~^» 
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C The calculation of density matrix elements for 
noncommutative theory 

The density matrix elements for the noncommutative theory are calculated in this ap- 
pendix. In the Boltzmann equation (13), the commutator [H int (0),V^] has the form 

[7^„ t (0),2^(k)] = J dqdc[d v d v '{27,f5\c( + v '-^-v){M e 1 + M d 2 + Ml + Ml) 

x[^(g')&r(gi)4(p')s(A;)2p /0 (27r) 3 ^5 3 (p - k) 

-^(g')^(g)ai(A:)a s (p)2p' (27r) 3 5, s ,5 3 (p' - k)], (90) 

where Ai® are the noncommutative parameter dependent parts of the amplitudes of the 
diagrams shown in Fig. 2. In this figure the box shows vertex in the noncommutative 
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q+p 




q+p 



q-p 





Figure 2: Diagrams for noncommutative Compton scattering 

space. These amplitudes can be calculated straightforwardly as follows 
M e =M e 1 + M 9 2 + M e 3 + M e 4 

= ^-uA<l') \UAp'W0(q + P ) + q'BesW) (rf + j> + m)^ s (p) 
op ■ q l\ / 

+is>{p'){i + i> + m) (jt s {p){q + p)9q + q9e s (p)p' ) u r (q) 

+ W~q UAq ' ] W s {p)q ' 6{q ~ P ' ] ~ q ' 6esW ) W " /+ m)UP ' ] 

+ts (P)d ~ /+ m) [W){q - p')9q - q9e s ,( P ')t )]u r (q), (91) 

where p9q = p^O^q^. Therefore one finds 

<([7L*(0),I>y(k)]) =^^(0) |rfq(27r) 3 ^(^p^(k)-^(k)) 

xv(?) [<zM*0 (Ws (k) + ^ s (&)#) + q6e s {k) (^(A;)# + u r (g), (92) 

After some calculations the time evolution of the density p^- generally is given by 
d , e 2 f , n e (q 



dq- 



(S is p s/j (k) - S js >p is (\x.)) 



dt rijy ' 16m J " k-q 
x (q6e s >(k)(16k ■ qq ■ e s ) + q9e s (k)(lQk ■ qq ■ e s >)\, 



(93) 



where for the special cases p\ 2 and p 21 one obtains 



Pi2(k) = y rfqn e (q)^2(g6»eig • e 1 + g6»e 2 g • e 2 )p 12 - (q9e 2 q ■ e l + q9e x q ■ e 2 )p u 



+ (q9e 2 q ■ t x + q9e x q ■ e 2 )p 22 ), 



(94) 



and 



p 2 i(k) = ~2^o J rf q^e(q)(2(g6»eig • e 1 + q9e 2 q ■ e 2 )p 21 - {q9e 2 q ■ e x + q9e x q ■ e 2 )p u 



+ 



(q9e 2 q ■ e 1 + q9e x q ■ e 2 )p 22 ), 



(95) 
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D The calculation of density matrix elements for 
magnetic field background 

In this appendix we calculate the density matrix elements in the presence of a back- 
ground magnetic field. Vanishing of the refractive term on the right hand side of (13) 
forces us to go to the next order and calculate the damping term of the Boltzmann 
equation (13). As it has been clarified in [3], using the expectation values of creation 
and annihilation operators calculated in this reference, the damping term can be written 
in terms of amplitudes in the following form 



Up to the linear order of B a p, the evaluating of the product of matrix elements yields 
Miq'r'jps'^qr, fcsi) x M(qr,ks' 2 ,q'r',ps 2 ) | B -term 



where the B-term shows the terms that are linear in terms of the background magnetic 
field. For relevant cosmological situations one can consider the cold plasma of electrons 
and photons that their kinetic energies are small compared to the electron mass i.e. 
p,k and q <^ m. As mentioned in [3], if the electron and photon temperatures are 
comparable which means p,k <C q, then various functions in (96) can be expanded in 
terms of p/q and q/m. Hence by computing the traces in (97), for the leading order, 




-(27r) 3 5 3 (0) J dqdpdcfM{q'r',ps' 1 ,qr,ks 1 )M{qr,ks 2 ,q'r',ps 2 ) 

x(27r) 4 5 4 (g' + p - q - fc)[n e (q)5 S2S ' i (5 isi p4 J (k) + 5 js > 2 p isi (\t)) 
-2n e {q)5 isi 5 js ' 2 p s ' iS2 {p)]. 



(96) 
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the matrix element is simplified as 



M(q'r',ps[,qr,ks 1 ) xM(qr, ks' 2) q'r' \ps 2 ) |b 

ie 4 ( 1 1 



-term 
1 



2m 2 \q.k q.p J \(q.k) 2 (q-p) 

q.e S2 (p)q.e si (fc)g.e si {p)q.e s > 2 (k) - q.e s > 2 (k)q.e S2 (p)q.e si (A;)g.e s ; (p) 
+q.e Sl (%.e s / (p)q.e S2 {p)q.e s > 2 (k) - q.e^ (p)q.e Sl (k)q.e S2 (p)q.e s > 2 (k) 
+0(k,p)], (98) 

where g M = —eB^ v q v '. Now by substituting (98) in (96) and integrating over the q', 
gives functions such as n e (q + (k — p)) and S(k —p + E(q) — E(q + k — p)) where can 
be expanded as follows 



ra e (q+ (k-p)) ~ n e (q) 



l _ (k-p).(q-mv) _ (k-p) 2 



2mT P 



(99) 



and 



S(k —p + E(q) -E(q + k- p)) ~ S(k -p)+ ( k P)-^S(k p) + 

m op 

Inserting the above equations into the (96) leads to the density matrix p i2 and p 2 i- 
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